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Abstract. Electric field dynamics at a positive ion imbedded in an electron 
gas is considered using a semiclassical description. The dependence of the field 
autocorrelation function on charge number is studied for strong ion-electron 
coupling via MD simulation. The qualitative features for larger charge numbers 
are a decreasing correlation time followed by an increasing anticorrelation. 
Stopping power and related transport coefficients determined by the time integral 
of this correlation function result from the competing effects of increasing initial 
correlations and decreasing dynamical correlations. An interpretation of the MD 
results is obtained from an effective single particle model showing good agreement 
with the simulation results. 
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1. Introduction 

The total electric field at a particle in a plasma determines the dominant radiative and 
transport properties of that particle. The properties of fields due to positive ions at a 
positive particle have been studied in some detail for both the static distribution 
of fields Q and the dynamics of the electric field autocorrelation function |j. 
The latter poses a real challenge since finite charge on the site at which the field is 
considered precludes the use of standard linear response methods. The corresponding 
study of fields at a positive ion due to electrons has been considered more recently 
for the simplest case of a single ion of charge number Z in an semiclassical electron 
gas. The static properties (electron charge density, electron microfield distribution) 
have been discussed in some detail elsewhere [||. Here, attention is focused on 
the dynamics via the electric field autocorrelation function. The case of opposite 
sign electron fields at a positive ion is qualitatively different from same sign ion 
fields, since in the former case electrons are attracted to the ion leading to strong 
electron-ion coupling for the enhanced close configurations. It is difficult a priori 
to predict the qualitative features of the correlation function due to this inherent 
strong coupling. Consequently, the analysis here has been based on MD simulation of 
the correlation function followed by an attempt to model and interpret the observed 
results. The simulations represent classical mechanics for Coulomb interactions with 
the ion-electron potential modified at short distances to represent quantum diffraction 
effects. The details have been discussed elsewhere Q and will not be repeated here. 
There are only three dimensionless parameters: the charge number of the ion, Z, the 
electron-electron coupling constant T, and the de Broglie wavelength relative to the 
interelectron distance, 5. The electron-ion coupling is measured by the maximum value 
of the magnitude of the regularized ion-electron potential at the origin, a = ZT/5. In 
this brief report results are reported for T = 0.1 and a = 0.25Z, with Z = 8,20,30, 
and 40. The corresponding density and temperature are n = 2.5 x 10 22 cm~ 3 and 
T = 7.9 x 10 5 K. 

The primary observations from the simulations of the field autocorrelation 
function for increasing charge number are: 1) an increase in the initial value, 2) a 
decrease in the correlation time, and 3) an increasing anticorrelation at longer times. 
The stopping power for the ion by the electron gas is proportional to the time integral 
of the field autocorrelation function in the limit of large ion mass This same 
integral determines the self-diffusion and friction coefficients in this same limit Q . The 
nonlinear dependence of these properties on Z is therefore the result of competition 
between the increase of the integral due to 1) and the decrease due to 2) and 3). 
The results from simulation show that the latter two dynamical effects dominate the 
former static effect. To interpret this, a simple model for the field correlation function 
is proposed such that the initial correlations are given exactly, but the dynamics is 
determined approximately from a single electron-ion trajectory. The model reproduces 
well the MD simulation results and suggests an interpretation of 2) and 3). 

2. MD simulation of field dynamics 

The system considered consists of N e electrons with charge — e, an infinitely massive 
positive ion with charge Ze placed at the origin, and a rigid uniform positive 
background for overall charge neutrality. The regularized electron-ion potential is 

chosen to be V(r) = -Ze 2 (l - e~ r / s ) /r where 8 = {2nh 2 /m e k B T) 1/2 is the 
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Figure 1. Field autocorrelation function for Z = 8,20,30,40 at T = 0.1 and 
5 = 0.4. 



electron thermal de Broglie wavelength. For values of r >> 5 the potential becomes 
Coulomb, while for r « 5 the Coulomb singularity is removed and V(r) — > —Ze 2 /5. 
This is the simplest phenomenological form representing the short range effects of 
the uncertainty principle |7j . Dimensionless variables are based on scaling coordinates 
with the average electron-electron distance, r , defined in terms of the electron density 
n e by 47rn e rQ/3 = 1, and scaling time with the inverse electron plasma frequency. The 
electron electric field at the ion is obtained from the total regularized potential 

N 

E = -V V({r l0 }) = Y,e(r t0 ) (1) 

2 = 1 

where r^o = i"i — ro is the position of the i th electron relative to the ion, and 

N ^ 

V({r i0 }) = ^2V(r i0 ), e (r i0 ) = e^(l-(l + ^)e-^/ 5 ) (2) 

1=1 l ® 

The dimensionless field autocorrelation function is denned by 

C(t) = 4<E(t)-E>. (3) 

The brackets denote an average over the classical Gibbs ensemble for the composite 
electron-ion system at equilibrium. 

The results for C(t) from MD simulation are shown in Figure 1 for Z = 8, 20, 30, 
and 40. The initial value increases approximately as a third order polynomial in Z g. 
The correlation time t c is defined to be the time at which the correlation function first 
goes to zero, C{t c ) — 0, and is seen to decrease as Z increases. Finally, in all cases 
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Figure 2. Stopping power S(v)/v for Z = 8,20,30,40 at V = 0.1 and <5 = 0.4. 



there is anticorrelation (C(i) < 0) for t > t c . The physical basis for the increase in 
the initial value is easy to understand. As Z increases, the electron density near the 
ion increases and the magnitude of the field increases for these more probable closer 
configurations. An explanation for the correlation time and anticorrelation is more 
difficult, and is the objective of the following sections. First, some consequences of 
this behavior are illustrated. 

3. Stopping power, friction, and self-diffusion 

The case of an infinitely massive ion considered here leads to exact relationships 
between transport coefficients characterizing three physically different phenomena: 
1) the low velocity stopping power S for a particle injected in the electron gas, 2) the 
friction coefficient £ for the resistence to a particle being pulled through the gas, and 
3) the self-diffusion coefficient D of a particle at equilibrium with the gas. The exact 
relationship is |6| 

/•OO 

m £ = (/3£>) _1 = S(v)/v = f3Z 2 r- 4 / dtC(t) (4) 

Jo 

This Green-Kubo representation allows determination of these transport properties 
from an equilibrium MD simulation. Previous simulations of stopping power have 
studied the nonequilibrium state of the injected particle, measuring directly the energy 
degradation ^ . As discussed below, (Q) provides the basis for additional interpretation 
of the results. 

Figure 2 shows the dimensionless stopping power obtained from the results of 
Figure 1 as a function of Z. Also shown is the Born approximation cx Z 2 , valid for 
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small Z. Previous simulations || and some experiments || suggest a crossover at 
larger Z to a weaker growth oc Z 2 . The data in Figure 2 has been fit to a such a 
crossover function showing consistency with these earlier results. This behavior is 
somewhat puzzling in light of the strong growth of the initial value C(0) « Z 3 at 
large Z. Making this explicit, the stopping power can be expressed as 

POO 

S(v)/v tx Z 5 / dtf(t), f(t) = C(t)/C(0) (5) 

Jo 

Evidently, the dynamical effects of the normalized correlation function f(t) decrease 
the stopping power as ~ Z~ 3 5 for large Z. A possible explanation for this is given in 
the following section. 



4. Single particle dynamics 

Consider again the initial covariance C (0) for which three representations can be given 
C (°) = J; / dre ^- ^(r)e(r) + J- J dr' fc (r,r')e(r') 
= db /^aWV-e(r) 

= ^ / dr 9ie(r)e mf (r) e (r) (6) 

The first equality expresses the covariance in terms of the one and two electron 
correlations with the ion, gi e (r) and <?i e (r, r') respectively. The second equality 
exploits the relationship for the field to the Gibbs factor (3U({rio}) = ZTV({rio}) 
and an integration by parts. This second representation requires only the one electron 
correlation function. The third representation is obtained from the second by an 
integration by parts to identify the mean force field e TO / (r) = V In g ie (r) . The third 
equality of (||) is similar to the first with apparent neglect of the two electron-ion 
correlations. However, these latter contributions are incorporated exactly in the mean 
force field e m / (r). This suggests a corresponding model for finite times 

C ( t ) ^ J dvdr<f>(v)g le (r)e m/ (r) e (r(t)) (7) 

where </>(i>) is the normalized Maxwellian and r(t) is a single electron trajectory in the 
presence of the ion, generated by the potential associated with 3i e (r) for stationarity. 
Further details and a justification based on the Vlasov equation will be given elsewhere. 
Clearly, the initial covariance is given exactly by this approximation. For practical 
purposes gi e {r) and the associated field and potential are represented here by the non- 
linear Debye-Huckel approximation with effective charge number and screening length 
adjusted to fit the MD data. 

Figure 3 shows a comparison of (0) with the MD results of Fig. 1 for Z = 8 and 
Z = 30. The agreement is very good, indicating that the dominant Z dependence is 
captured by the initial correlations and the single particle dynamics. 



5. Discussion 



The good agreement of the simple model in the previous section suggests that the 
decrease in correlation time and build up of anticorrelation can be understood in 



Electron dynamics at a positive ion 



6 



A 



15- 




• O MD 


10- 




Model 


5- 


\ Z=8 




0- 








' I ' I ' I ' I 



0.0 0.5 1.0 1.5 2.0 



CO t 
P 



Figure 3. Comparison of C(t) from MD and from EqjT] for Z = 8 and 30 at 
T = 0.1 and 8 = 0.4. 

terms of one electron dynamics. As Z increases the probability of close electron 
ion configurations increase. The electron is subjected to greater acceleration toward 
the ion and the time to reverse its sign decreases. This is the effect of decreasing 
correlation time. Once the field has reversed sign there is anticorrelation. Since the 
closer configurations imply larger field values the magnitude of this anticorrelation also 
increases with increasing Z . It remains to quantify this picture but the single electron 
dynamics appear to provide qualitative confirmation. A more complete discussion will 
be provided elsewhere. 
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